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Elastic objects across a wide range of scales deform under local changes of their intrinsic properties, 
yet the shapes are glocal, set by a complicated balance between local properties and global geometric 
constraints. Here, we explore this interplay during the inversion process of the green alga Volvox, 
whose embryos must turn themselves inside out to complete their development. This process has 
recently been shown [S. Hohn et al, Phys. Rev. Lett. 114, 178101, (2015)] to be well described by 
the deformations of an elastic shell under local variations of its intrinsic curvatures and stretches, 
although the detailed mechanics of the process have remained unclear. Through a combination of 
asymptotic analysis and numerical studies of the bifurcation behavior, we illustrate how appropriate 
local deformations can overcome global constraints to initiate inversion. 


I. INTRODUCTION 


The shape of many a deformable object arises through 
the competition of multiple constraints on the object: 
this competition may be between different global con¬ 
straints, such as in Helfrich’s analysis [1] of the shape of 
a red blood cell (where intrinsic curvature effects coexist 
with constrained membrane area and enclosed volume). 
It may also be the competition between local and global 
constraints. Such deformations, which we shall term 
glocal, arise for example in origami patterns [2] (where 
local folds must be compatible with the global geome¬ 
try). They are of considerable interest in the design of 
programmable materials [3] at macro- and microscales, 
where one asks: can a sequence of local deformations 
overcome global constraints and direct the global defor¬ 
mations of an object? 

This is a problem that, at the close of their devel¬ 
opment, the embryos of the green alga Volvox [4] are 
faced with in the ponds of this world. Volvox (Fig. la) 
is a multicellular green alga belonging to a lineage (the 
Volvocales) which has been recognized since the time of 
Weismann [5] as a model organism for the evolution of 
multicellularity, and which more recently has emerged as 
the same for biological fluid dynamics [6]. The Volvocales 
span from unicellular Chlamydomonas, through organ¬ 
isms such Gonium, consisting of 8 or 16 Chlamydomonas- 
like cells in a quasi-planar arrangement, to spheroidal 
species [Pandorina and Pleodorina) with scores or hun¬ 
dreds of cells at the surface of a transparent extracellular 
matrix (ECM). The largest members of the Volvocales 
are the species of Volvox, which display germ-soma dif¬ 
ferentiation, having sterile somatic cells at the surface of 
the ECM and a small number of germ cells in the interior 
which develop to become the daughter colonies. 

Following a period of substantial growth, the germ cells 
of Volvox undergo repeated rounds of cell division, at the 
end of which each embryo (Fig. lb,e) consists of a few 
thousand cells arrayed to form a thin spherical sheet [4]. 
These cells are connected to each other by the remnants 
of incomplete cell division, thin membrane tubes called 
cytoplasmic bridges [7, 8]. The ends of the cells whence 


emanate the flagella, however, point into the sphere at 
this stage, and so the ability to swim is only acquired 
once the alga turns itself inside out through an opening 
at the top of the cell sheet, called the phialopore [9-11]. 

Of particular interest in the present context is the cru¬ 
cial first step of this process, the formation of a circu- 



FIG. 1. (color online). Vbfooa; invagination and elastic model. 

(a) Adult Volvox, with somatic cells and one embryo labelled. 

(b) Volvox embryo at the start of inversion, (c) Mushroom¬ 
shaped invaginated Volvox embryo, (d) Cell shape changes to 
wedge shapes and motion of cytoplasmic bridges (CB) bend 
the cell sheet. Red line indicates position of cytoplasmic 
bridges. (e,f) Cross-sections of the stages shown in panels 
(b,c). Cell shape changes as in (d) occur in the marked re¬ 
gions. (g) Geometry of undeformed spherical shell of radius 
R and thickness h. (h) Geometry of deformed shell. Scale 
bars: (a) 50 pm, (e,f) 20 pm. False color images obtained from 
light-sheet microscopy provided by Stephanie Hohn and Au¬ 
relia R. Honerkamp-Smith. 
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lar invagination in so-called ‘type B’ inversion (Fig. lc,f) 
followed by the engulfing of the posterior by the anterior 
hemisphere [11, 12]. (This scenario is distinct from ‘type 
A’ inversion in which the initial steps involve four lips 
which peel back from a cross-shaped phialopore.) The 
invaginations of cell sheets found in type B inversion are 
very generic deformations during morphogenetic events 
such as gastrulation and neurulation [13-16], but, in an¬ 
imal model organisms, they often arise from an intricate 
interplay of cell division, intercalation, migration and cell 
shape changes. Modelling these therefore requires cell- 
based models, as pioneered by Odell et al. [17], but 
simpler models of simpler morphogenetic processes are 
required to elucidate the underlying mechanics of these 
problems [18]. Inversion in Volvox is, however, driven 
by active cell shape changes alone: inversion starts when 
cells close to the equator of the shell elongate and be¬ 
come wedge-shaped [12]. Simultaneously, the cytoplas¬ 
mic bridges migrate to the wedge ends of the cells, thus 
splaying the cells locally and causing the cell sheet to 
bend [12] (Fig. Id). Additional cell shape changes have 
been implicated in the relative contraction of one hemi¬ 
sphere with respect to the other in order to facilitate 
invagination [19]. After invagination, the bend region ex¬ 
pands, allowing the posterior hemisphere to invert fully. 

At a more physical level, it has been shown re¬ 
cently that the inversion process is simple enough to be 
amenable to a mathematical description [19]: the defor¬ 
mations of the alga are well reproduced by a simple elas¬ 
tic model in which the cell shape changes and motion of 
cytoplasmic bridges impart local variations of intrinsic 
curvature and stretches to an elastic shell [19]. The asso¬ 
ciated mechanics have remained unclear, however. Here, 
we perform an asymptotic analysis at small deformations 
to clarify the geometric distinction between deformations 
resulting from intrinsic bending and intrinsic stretching, 
respectively. A numerical study of the bifurcation be¬ 
havior further serves to illustrate how a sequence of local 
deformations can achieve invagination, and how contrac¬ 
tion complements bending in this picture. 


II. ELASTIC MODEL 

Following Hohn et al. [19], we inscribe Volvox inver¬ 
sion into the very general framework of the axisymmet- 
ric deformations of a thin elastic spherical shell of radius 
R and thickness h R under variations of its intrin¬ 
sic curvature and stretches. The undeformed, spherical, 
configuration of the shell is characterized by arclength s 
and the distance of the shell from its axis of revolution, 
ro(s) (Fig. Ig). To these correspond arclength S{s) and 
distance from the axis of revolution r(s) in the deformed 
configuration (Fig. Ih). The undeformed and deformed 
configurations are related by the meridional and circum¬ 
ferential stretches, 

Ms) = ^ and /^(s) = . (1) 

ds ro(s) 


(These definitions do not require that the undeformed 
configuration be spherical, and apply for the deforma¬ 
tions of any axisymmetric object.) These define the 
strains 

Es = fs- /°, E^ = U- /o, (2) 

and curvature strains 

Ks = fsKs - /°K°, - f^K^, (3) 

where Ks and denote the meridional and circumfer¬ 
ential curvatures of the deformed shell. The intrinsic 
curvatures and stretches introduced by /°, and k®, 
extend Helfrich’s work on membranes [1]. The deformed 
configuration of the shell minimises an energy of the 
Hookean form [20-22] 

^ I + 21 .A, ds 

+ i ds. (4) 

with material parameters the elastic modulus E and Pois¬ 
son’s ratio h'. In computations, we take h/R = 0.15 and 
v = 1/2 appropriate for Volvox inversion [19]. 

In general, deformations of the shell arise from a com¬ 
plex interplay of intrinsic stretches and curvatures, and 
the global geometry of the shell. To clarify these, we be¬ 
gin by considering two simple kinds of deformations, in 
which the competition is between two effects only. How 
these effects conspire in general we shall explore in the 
main body of the paper. 

1. Simple Deformations: Stretching and Bending 

The simplest intrinsic deformation is one of uniform 
stretching or contraction, which does not affect the 
global, spherical geometry of the shell. This corresponds 
to fg = f/p = f and Kg=Kp = 1/fR. With these intrin¬ 
sic stretches and curvatures, the original sphere deforms 
to a sphere of radius R'. Then fg = fp = R'/R, and so 
the strains are Eg = Ep, = R'/R—f. However, Kg = = 

1/R'. Thus fgKg = fcjyKcj, = = 1/^7 and so 

Kg = Kp = 0. The energy density is therefore propor¬ 
tional to {R'/R — M and is minimized for R' = fR, at 
which point £ = 0 (Fig. 2a). (Indeed, uniform contrac¬ 
tion is a homothetic transformation: the angles between 
material points are unchanged, and so there is no bending 
involved. In other words, the shell is blind to its intrinsic 
curvature on this spherical solution branch.) 

The intrinsic stretches and curvatures need not be com¬ 
patible in this way, however: suppose that /g = /° = /, 
but Kg = = tj j'R with f ^ f. The energy still 

has spherical minima of radius El! = fR, but now with 
£ fz 0 (Fig. 2a). This illustrates that, conversely, even 
if the equilibrium shape is spherical, the intrinsic curva¬ 
tures and stretches cannot straightforwardly be inferred 
from the resulting shape. 
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FIG. 2. (color online). Simple intrinsic deformations, (a) A 
sphere can be shrunk to smaller spheres of equal radii by both 
compatible and incompatible intrinsic deformations, (b) Con¬ 
traction of a circular region of radius 7? in a plane elastic sheet 
by a factor /. The boundary of this region is contracted to 
s = FR. (c) Numerical result for F (-I-) agrees with analytical 
calculation (9) (solid line). 


2. Simple Deformations: Stretching and Geometry 

To illustrate how the global geometry affects these de¬ 
formations, we consider contraction of a plane elastic 
sheet, with /° = /^ = / < 1 for s < i? (Fig. 2b). This 
does not involve any bending of the sheet, and, upon 
non-dimensionalising lengths with i?, the shell minimises 

poo 

s{ [r'(s) - /(s)] ^ -b [r(s)/s - /(s)] ^ 

-b 2^[r'(s) - /(s)] [r(s)/s - /(s)] | ds, (5) 

where 

= ( 6 ) 

The resulting Euler-Lagrange equation is 

This is a homogeneous equation, and the solution satis¬ 
fying the geometric conditions r(0) = 0 and r(s) ~ s as 
s —> 00 as well as continuity of r at s = 1 is 

Fs if s < 1 


s 

The constant F = r(l) is determined by the jump con¬ 
dition at s = 1, or, physically, by requiring the stress to 
be continuous across s = 1. This finally yields 

F=^[(l-^) + (l + ^)/]. (9) 

This simplified problem serves as a test case for numer¬ 
ical solution of the more general Euler-Lagrange equa¬ 
tions associated with (4). These boundary-value prob¬ 
lems can be solved numerically with the solver bvp4c of 



Matlab (The Mathworks, Inc.); our numerical setup of 
the governing equations otherwise mimicks that of [22]. 
In this particular example, the linear relationship in (9) 
is indeed confirmed numerically (Fig. 2c). Notice that 
the governing equation (7) is independent of the forcing 
applied away from s = 1; the solution is determined by 
geometric boundary conditions. 


III. RESULTS 


The most drastic cell shape changes at the start of in¬ 
version occur when cells in a narrow region close to the 
equator become wedge-shaped (Fig. Id). These are ac¬ 
companied by motion of the cytoplasmic bridges to the 
thin tips of the cells to splay the cell sheet and drive its 
inward bending. For this reason, Hohn et al. [19] started 
by considering a piecewise constant functional form for 
the intrinsic curvature, in which this curvature took neg¬ 
ative values in a narrow region close to the equator. It 
was found, however, that with this ingredient alone the 
energy minimizers could not reproduce the mushroom 
shapes adopted by the embryos in the early of stages of 
inversion (Fig. lc,f), producing instead a shape cinched in 
at those points - the so-called ‘purse-string’ effect. How¬ 
ever, analysis of thin sections had previously revealed 
that the cells in the posterior hemisphere become thin¬ 
ner at the start of inversion [12]. When the resulting 
contraction of the posterior hemisphere was incorporated 
into the model, it could indeed reproduce, quantitatively, 
the shapes of invaginating Volvox embryos. 

Hohn et al. have thus identified two different types of 
active deformations that contribute to the shapes of in¬ 
verting Volvox at the invagination stage: first, a localized 
region of active inward bending (corresponding to nega¬ 
tive intrinsic curvature), and second, relative contraction 
of one hemisphere with respect to the other. We shall 
focus on these two types of deformation in what follows 
and clarify the ensuing elastic and geometrical balances. 


A. Asymptotic Analysis 


We start by seeking equilibrium configurations in the 
limit of a thin shell, h <C i?. In this limit, the shapes 
(Fig. 3a,b) corresponding to contraction or (pure) invagi¬ 
nation (by which we mean, here, deformations driven by 
a region of high intrinsic curvature only) result from the 
matching of spherical shells of different radii or disparate 
relative positions (Fig. 3c,d). Deviations from these outer 
solutions are localized to an asymptotic inner layer of 
non-dimensional width S about x = X, where x = s/R is 
the angle that the normal to the undeformed shell makes 
with the vertical (Fig. 3e). Here, we consider an incipient 
deformation where the normal angle f3{x) to the deformed 
shell deviates but slightly from its value in the spherical 
configuration, viz I3{x) = x + b{x), with 6 <C 1. 
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FIG. 3. (color online). Asymptotic analysis of invagination and contraction, (a) Numerical shape resulting from contracting 
the posterior to a radius Vp < R. (b) Numerical “hourglass” shape resulting from pure invagination, (c) Geometry of 
contraction with posterior radius rp < R, resulting in upward motion of the posterior by a distance d. (d) Geometry of pure 
invagination solution, (e) Asymptotic geometry: in the limit h R, deformations are localized to an asymptotic inner layer 
of width 5 about x = X, where x = s/R is the angle that the undeformed normal makes with the vertical. In the deformed 
configuration, this angle has changed to P{x). (f) Asymptotic invagination: upward motion of posterior by a distance d requires 
inward deformations scaling as {Sd)^^^ in the inner layer of width 6. (g) Relation between preferred curvature k and width of 
invagination A for a given amount of upward posterior motion d, from asymptotic calculations, (h) Inward rotation of midpoint 
of invagination with, and without contraction, from asymptotic calculations. 


1. Geometric Considerations 

We begin by clarifying the geometric distinction be¬ 
tween contraction and invagination. The radial and ver¬ 
tical displacements obey 

u'r = fs cos/3 — cosx = —6sin AT -|- 0(^Sb, 6^), (10a) 
u'^ = fs sin/3 — sin a; = bcosX + 0{5b, 6^), (10b) 

where dashes denote differentiation with respect to cc, 
and where we have assumed the scaling fg = l + 0{Sb) 
which we shall derive presently. Let d denote the (non- 
dimensional) distance by which the posterior moves up. 
Matching to the outer solutions requires the net displace¬ 
ments Ur and C4, obtained by integrating (10) across the 
inner layer, to obey 

= dsinX, = —dcosX, (Ha) 

C/O) = 0, up = -d, (lib) 

where the superscripts (c) and (i) refer, respectively, to 
the solutions corresponding to contraction and (pure) in¬ 
vagination. In the case of contraction, (10) and (11a) give 
the scaling b^^'> ^ d/S. If there is no contraction, however, 
(10) and (11b) imply that the leading-order solution does 
not yield any upward motion of the posterior, which is as¬ 
sociated with a higher-order solution only. This suggests 
that the appropriate scaling is ~ {d/Sy^^, which we 
shall verify presently. 

Our assumption 6 1 thus translates to d 6 . 

Hence, in the invagination case, upward motion of the 
posterior requires comparatively large inward displace¬ 
ments of order ((5d)^/^ 3> d (Fig. 3f). This asymptotic 


difference of the deformations corresponding to contrac¬ 
tion and invagination arises purely from geometric ef¬ 
fects; it is the origin of the ‘purse-string’ shapes found 
by Hohn et al. in the absence of contraction [19]. 


2. Elasto-Geometric Considerations 


Here, we discuss the detailed solution for pure invagi¬ 
nation. Upon non-dimensionalising distances with R and 
stresses with Eh, the Euler-Lagrange equations of (4), 
derived in the appendix, can be cast into the form 


fs X sin X tan /3 — cos /3 ^1 — vjd' — sin /3 cosec xj 

— ^/3'sinx-I-Jz(sin/3 — sina:)^ = k°{x), (12a) 

fusec^dsina;') — A — izU = 0, (12b) 

da; V / 

with the small parameter 


£ 


2 


1 

12 ( 1 - 1 / 2 ))^ 


< 1 . 


(13) 


In these equations, E is the non-dimensional meridional 
stress, and A = is the dimensionless hoop strain. The 
contribution from the intrinsic curvature is 


k°{x) = 


|i/Ks(a) cos/3 



(14) 


The equations are closed by the geometric relation 


d_ 

dx 


{^A sin Pj = fs cos /3 


cos a;. 


(15) 
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Introducing 7 = scaling gives the leading bal¬ 

ances U ~ £^ 7 /( 5 ^, I^/S ~ A, and A/5 ~ 7 in (12,15). 
Hence <5 ~ £^/^, and we define an inner coordinate ^ via 
X = X + 5^. We also introduce the expansions 

/3 = X-I- 7(60 + 7 ^ 1 -b 7^52 H-), (16a) 

A = 5^(ao + ^ai-\ -), (16b) 

X = ( 5^7 cot X ((To + TO"! -!-•••). (16c) 

This further proves the scaling fs = l + 0{Sj) that we 
have assumed previously. 

The pure invagination configuration is forced by in¬ 
trinsic curvature that differs from the curvature of the 
undeformed sphere in a region of width A about x = X, 
where = —k. Writing A = we thus have, at 

leading order, 


At next order, (11b) is a system of two linear equations 
for two integrals, with solution 

/ OO POO 

60 d^ = 2 sin X, / 6 ^ d^ = — cos X. (23) 

-00 J —OO 

We note in particular that the resulting condition on the 
leading-order solution has only arisen in the second-order 
expansion of the matching conditions. Similarly, at order 
0 (( 57 ^), we find 


bohi d^ = 0. 


The leading-order problem is thus 




where k = and where H denotes the Heavi¬ 

side function. Thus 

fcO(^) = sinX, (18) 

where = X [d -b lyl) — 5 — Id)]. We note 

that 7 ^ ^ 'yS provided that e (which we shall assume 
to be the case); thus we may set a; = X to the order at 
which we are working. Expanding (12,15), we then find 


(To - = X°( 5 ), (To-ao = 0 , 


Jn — —bn 


at lowest order, where dashes now denote differentiation 
with respect to At next order, 

(Ti -b (T 060 secX cosec X — 6 " = 0, (20a) 

(t( -b ^ (^o^o) sin X sec X — tti = 0, (20b) 

with a'l = —bx — \bn cotX. We are left to determine the 
matching conditions by expanding ( 10 ) to find 

u/ = —S'fbo sinX — 5^"^ {pi sinX -b cosX) 

— ( 57^(62 sinX -b hnbi cosX — IfepSinX), (21a) 
u/ = S'-fbo cos X -b ^ 7 ^ (bi cos X — I^q sin X) 

-b 5j^(b2 cosX — bobi sinX — Ifep cosX), ( 21 b) 

up to corrections of order Applying (11b), 

at lowest order, we find 


bo d^ = 0. 


with matching conditions (22) and the first of (23). Sym¬ 
metry ensures that the first of (22) is satisfied. After a 
considerable amount of algebra, the first of (23) reduces 
to a relation between K and A, 


8-\/2 sin X 


1 _|_ e-^/'^(^(\/2A - 1 ) sin ^ - cos 


■ (26) 


This function exhibits a global minimum at d « 2.44 
(Fig. 3g). This is a first indication that narrow invagina¬ 
tions are more efficient than those resulting from wider 
regions of high intrinsic curvature, a statement that we 
shall make more precise later. 

Symmetry also implies that there is no inward rotation 
of the midpoint of the invagination at this order. Rather, 
inward folding is a second-order effect, for which we need 
to consider the second-order problem, 

b'i” + = I ^ (^o^o) -\bi^ cot X, (27) 

with matching conditions (24) and the second of (23). 
The rotation of the midpoint of the invagination is thus 


A/ 36 ) = (H 6 )(d) cos. 


where i? 6 )(yl) is determined by the solution of (27). The 
detailed solution reveals that 


h 6 )(A) 


2i/2e 4e sin ^ — e sin — 3 sin -b ^e — 1^ cos 
5 e73 -b (v^d — 1 ) sin ^ — cos ^ 
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FIG. 4. (color online). Setup for numerical calculations, fol¬ 
lowing [19]. (a) Geometrical setup: the intrinsic curvature k® 
of a spherical shell of undeformed radius R differs from the un¬ 
deformed curvature in the range Amax > s > Amax — A, where 
s is arclength. Posterior contraction is taken into account by 
a reduced posterior radius Vp < R. These intrinsic curvature 
and contraction result in deformations that move up the pos¬ 
terior pole by a distance d. (b) Corresponding functional form 
of in the bend region, = —k < 0. (c) Form of the intrin¬ 
sic stretches fs,f^ for posterior contraction, (d) Functional 
form of for posterior contraction. 


but the geometric factor in (28) is the main point: this 
factor resulting from the global geometry of the shell 
hampers the inward rotation of the midpoint of the in¬ 
vagination. (This is as expected: by symmetry, invagina¬ 
tion at the equator, where cosX = 0, yields no rotation.) 

An analogous, though considerably more straightfor¬ 
ward, calculation can be carried out for contraction: non- 
dimensionally, upward posterior motion by d requires 
= 1 — d ioT X < X, and leads to 




1 d 
2^2 £ 1 / 2 ' 


(30) 


At this order, the above solutions for pure invagination 
and contraction can be superposed; in particular, the so¬ 
lutions at order 0 ( 7 ^) have the same symmetry, and so 
(24) is satisfied. For contraction, there is thus no geo¬ 
metric obstacle to inward folding (Fig. 3h). Contraction 
is thus not only a means of creating the disparity in the 
radii of the anterior and posterior hemispheres required 
to fit the partly inverted latter into the former, but also 
drives the inward folding of the invagination, by breaking 
its symmetry. In Volvox inversion, this symmetry break¬ 
ing is at the origin of the formation of the second passive 
bend region highlighted by Hohn et al. [19] to stress the 
non-local character of these deformations. 


B. Bifurcation Behaviour 

The asymptotic analysis has shown that the coupling of 
elasticity and geometry constrains small invagination-like 



FIG. 5. (color online). Bifurcation behaviour of invagination 
solutions. Solution space for Amax = I.IR and Vp = R-. each 
line shows the relation between k and d at some constant 
A. A critical branch (at A = A*) separates different types of 
branches. Branches with A > A* feature two extrema; the 
resulting spinodal curve (thick dashed line) defines a critical 
point. Insets illustrate representative solution shapes. See 
text for further explanation. 


deformations both locally and globally, but that contrac¬ 
tion can help overcome these global constraints. These 
ideas carry over to larger deformations of the shell, which 
must however be studied numerically. For this purpose, 
we extend the setup of [19], motivated by direct obser¬ 
vation of thin sections of fixed embryos: the intrinsic 
curvature differs from that of undeformed sphere in 
the range Amax > s > Amax — A of arclength along the 
shell (Fig. 4a). In this region of length A, = —k, where 
k > 0 (Fig. 4b). This imposed intrinsic curvature results 
in upward motion of the posterior pole by a distance d. 


1. Stability Statements 

Our first observation is that, at fixed Amax, more than 
one solution may arise for the same input parameters 
(A:, A). Further understanding is gained by considering, 
at fixed Amax and for different values of A, the relation be¬ 
tween k and d. The typical behaviour of these branches 
is plotted in Fig. 5. (The shapes eventually self-intersect; 
accordingly, these branches end but we expect them to be 
joined up smoothly to configurations with opposite sides 
of the shell in contact. The study of such contact config¬ 
urations typically requires some simplifying assumptions 
to be made [22], but we do not pursue this further, here.) 

At the distinguished value A = A*, a critical branch 
arises (Fig. 5). It separates two types of branches: first, 
those with A < A*, on which d varies mononotonically 
with k, and second, those with A > A*, where the re¬ 
lation between d and k is more complicated. At large 
values of A, these branches may have a rather involved 
topology involving loops. At values of A just above A*, 
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FIG. 6. (color online). Stability of invagination solutions, (a) Possible topologies of a double fold in the distinguished (fc, x) 
bifurcation diagram. Dashed branches are those that the results of [23] imply to be unstable, (b) For A > A,, S-shaped folds 
arise in the {k, x) diagram. From general theory [23], the middle part of the branch is unstable, while the outer parts are stable. 
An additional region of metastability is identified by the Maxwell construction, (c) Resulting picture: a region of unstable and 
metastable solutions expands underneath the critical point. 


however, there is a range of values of k for which there 
exist three configurations (Fig. 5). We note that the two 
outer configurations have dd/dk > 0, while the middle 
one has dd/dk < 0. The latter behaviour prefigures in¬ 
stability, which we shall discuss in more detail below. 
There are thus two points on these branches where fc, 
viewed locally as a function of d, reaches an extremum. 
The curve joining up these extrema for different values of 
A we shall term the ‘spinodal curve’. This curve, in turn, 
has a maximum at a point on the critical branch, which 
we shall call the ‘critical point’ and which is characterised 
by A* and the critical curvature, fc*. 

The stability of the configurations in Fig. 5 can be as¬ 
sessed by means of general results of bifurcation theory 
[23], used recently to discuss the stability of the buckled 
equilibrium shapes of a pressurised elastic spherical shell 
[22, 24]. If we let >f = —dEjdk denote the conjugate vari¬ 
able to k, the key result of [23] is that stability, at fixed A, 
of extremizers of the energy £ can be assessed from the 
folds in the (fc, >£■) bifurcation diagram. In particular, sta¬ 
bility can only change at folds in the bifurcation diagram. 
Expanding the bending part of the energy functional (4) 
for /° = /^ = 1 and = 0, we find 

•^max y \. 

ro ifsKs + i'Uk^ + k \ ds (31) 

‘max ^ 

with G = 7rE/i^/6(l — (The last term in the inte¬ 
grand is independent of the solution, and may therefore 
be ignored in what follows.) 

The two folds that arise in the (fc, d) diagram for 
A > A* (Fig. 5) are compatible a priori with four fold 
topologies in the (fc, >c) diagram (Fig. 6a). However, since 
a single solution exists for small k (at fixed A), the lowest 
branch must be stable. Further, since the branches do 
not self-intersect in the {k,d) diagram, they cannot self¬ 
intersect in the (fc, k) diagram either. The results of [23] 
imply that only the first topology in Fig. 6a is compatible 
with this, and so the fold is S-shaped and traversed up¬ 
wards in the (fc, >c) diagram. It follows in particular that 
the middle branch, with dd/dk < 0 is unstable, and that 


right branch is stable. (Numerically, one confirms that 
the branches are indeed S-shaped.) Thus the stability of 
the branches in this simple bifurcation diagram could also 
be inferred from the (fc, d) diagram (though, in general 
problems, as discussed in [23], different bifurcation dia¬ 
grams may suggest contradictory stability results). How¬ 
ever, the Maxwell construction of equal areas [25] can 
be applied to the {k, >c) diagram (Fig. 6b) to identify 
metastable solutions beyond the unstable branch. These 
stability considerations may appear rather technical, but 
they are in fact very natural: under reflection, the (d, k) 
diagram maps to the diagram of isotherms of a classical 
van-der-Waals gas, for which the middle branch is well 
known to be unstable [25]. Under this analogy, £ corre¬ 
sponds to the Gibbs free energy of the gas. 

This analysis cannot immediately be extended to the 
more exotic topologies that arise for A close to Amax 
(Fig. 5). We note however that part of these branches 
must be unstable, too: as above, a single solution exists 
for small d, and so the corresponding branch must be sta¬ 
ble. The first fold must be traversed upwards, and the 
first branch with dd/dk < 0 is thus unstable, as above. 

An analogous analysis can be carried out for defor¬ 
mations that vary A while keeping k fixed: for k > k^^ 
the (d. A) diagram is monotonic, but this ceases to be 
the case for k < k^,. As above, the stability can be in¬ 
ferred from the (A, d) diagram, and the middle branch 
with 9d/i9A < 0 is unstable, too. 

The picture that emerges from this discussion is the 
following: solutions in a region of parameter space under¬ 
neath the critical point bounded by the spinodal curve 
are unstable; a band of solutions on either side of this re¬ 
gion and below the critical point are metastable (Fig. 6c), 
both to perturbations varying k and to perturbations 
varying A. If invagination, driven by a localized region 
of active bending, is to be stable, it must move around 
the critical point: if it were to enter the unstable region, 
the shell would flip back and forth between the ‘shal¬ 
low’ and ‘deep’ invagination states on either side of the 
unstable region and suffer large strains in the process. 

















FIG. 7. (color online). Contraction and the critical point. 
Trajectories of critical point in parameter space as Amax is 
varied, for different values of Cp. Thin dotted lines are curves 
of constant Amax. At constant Amax, increased contraction 
leads to decreased fc* and increased A*. 

(This makes this kind of instability different from the 
classical buckling instability of a rod or a ‘popper’ toy 
[26]: the latter is directed in that, once the instability 
threshold is crossed, the system will snap to the new pre¬ 
ferred configuration and remain there.) The need for a 
sequence of stable deformations to move around the crit¬ 
ical point rationalises the timecourse of invagination in 
Volvox: initially, a narrow band of cells undergoes cell 
shape changes, thereby acquiring a high intrinsic curva¬ 
ture. This region of cells then widens, moving around the 
critical point, whereupon the preferred curvature relaxes 
and posterior inversion can complete. 

2. Contraction and Criticality 

For different values of Amax, the critical point traces 
out a trajectory in parameter space, characterised by fc* 
and A* (Fig. 7). As Amax increases, /c* increases, while A* 
decreases. Thus the closer to the equator, the more dif¬ 
ficult invagination is, not only because there is less room 
to fit the posterior into the anterior, but also because 
a stable invagination requires narrower and narrower in¬ 
vaginations of higher and higher intrinsic curvature. 

We are left to explore how contraction affects the po¬ 
sition of the critical point, and hence the invagination. 
We introduce a reduced posterior radius rp < i? as 
in [19] (Fig. 4a), and modify the intrinsic curvatures and 
stretches accordingly (Fig. 4b,c,d). Numerically, we ob¬ 
serve that, at constant Amax, increasing contraction (that 
is, reducing Cp) decreases the critical curvature A:*, and 
increases A* (Fig. 7). Hence contraction aids invagination 
not only geometrially, but also mechanically: first, it al¬ 
lows invagination close to the equator (which would oth¬ 
erwise be prevented by different parts of the shell touch- 



FIG. 8. (color online). Examples of “preferred” deformation 
modes, which are solutions of (33). Of the modes shown, the 
middle one has the lowest elastic energy. 

ing), and second, it makes stable invagination easier, by 
reducing A:*. Thus, again, contraction appears as a me¬ 
chanical means to overcome global geometric constraints. 

3. Asymptotic Analogy 

In the asymptotic analysis in the previous section, we 
restricted ourselves to small deviations of the normal an¬ 
gle from the spherical configuration so that the prob¬ 
lem remained analytically tractable. While the lead¬ 
ing scaling balances remain the same for large rotations, 
the resulting non-linear “deep-shell equations” cannot be 
rescaled so that the dependance on X drops out [21]. 
Some further insight can, however, be gained in the 
shallow-shell limit AT ^ 1: in terms of the inner coor¬ 
dinate we write 

m=XB{0 and N(0=e5(0. (32) 

In the absence of forcing by intrinsic curvature or 
stretches, the leading-order balance is 

2S” = 1 - and B” = SB, (33) 

where dashes denote, as before, differentiation with re¬ 
spect to 

This balance arises also in the study of a spherical shell 
pushed by a plane [21]: at large indentations, the shell 
dimples and the plane remains in contact with it only in a 
circular transition region joining up the undeformed shell 
to the isometric dimple. With the matching conditions 
B —>• ±1 as ^ ±oo, (33) describe the leading-order 
shape of this transition region [21]. Remarkably, this 
deformation is independent of the contact force, which 
only arises at the next order in the expansion [21]. 

The appropriate boundary conditions for the invagi¬ 
nation case are B —> 1 as ^ ±oo, and non-constant 
solutions of (33) can indeed be found numerically (some 
solutions are shown in Fig. 8). In these modes, the defor¬ 
mations are, in a sense, large compared to intrinsic cur¬ 
vature imposed, making them geometrically ‘preferred’. 
Their existence lies at the heart of the bifurcation be¬ 
haviour discussed above. 


IV. CONCLUSION 

In this paper, we have explored perhaps the simplest 
intrinsic deformations of a spherical shell: elastic and 




























9 


geometric effects conspire to constrain deformations re¬ 
sulting from a localized region of intrinsic bending. Con¬ 
traction, a somewhat more global deformation, alleviates 
these constraints and thereby facilitates the stable transi¬ 
tion from one configuration of the shell to another. This 
rich mechanical behaviour makes a mathematically in¬ 
teresting problem in its own right, yet this analysis has 
implications for Volvox inversion and wider material de¬ 
sign problems. 

Experimental studies of Volvox inversion [12, 19] had 
revealed the existence of posterior contraction, and in¬ 
deed, the simple elastic model that underlies this paper 
can only reproduce in vivo shapes once posterior contrac¬ 
tion is included [19]. Of course, contraction is an obvious 
means of creating a disparity in the anterior and poste¬ 
rior radii required ultimately to fit one hemisphere into 
the other, but the present analysis reveals that, beyond 
this geometric effect, there is another, more mechanical 
side to the coin: if contraction is present, lower intrinsic 
curvatures, i.e. less drastic cell shape changes, are re¬ 
quired to stably invert the posterior hemisphere. This as¬ 
cribes a previously unrecognized additional role to these 
secondary cell shape changes (i.e. those occurring away 
from the main bend region): just as the shape of the 
deformed shell arises from a glocal competition between 
elastic and geometric effects, a combination of local and 
more global intrinsic properties allow inversion to pro¬ 
ceed stably. Thus, as we have pointed out previously, 
this mechanical analysis rationalises the timecouse of the 
observed cell shape changes, thereby lending further sup¬ 
port to the observation of Hohn et al. [19], that it is a 
spatio-temporally well regulated sequences of cell shape 
changes that drives inversion. Thus, the remarkable pro¬ 
cess of Volvox inversion is mechanically more subtle than 
it may initially appear to be. 

Intrinsic deformations that allow transitions of an elas¬ 
tic object from one configuration to another are of inher¬ 
ent interest in the material design context, and divide 
into two classes: first, snapping transitions for fast tran¬ 
sitions between states, studied in [3], and second, stable 
sequences of intrinsic deformations. The glocal behaviour 
of the latter is illustrated by the present analysis: in par¬ 
ticular, additional transformations such as contraction 
can increase the number of stable parameter paths be¬ 
tween configurations of the elastic object. In this mate¬ 
rial design context, non-axisymmetric deformations such 
as polygonal folds or wrinkles [27] could also become im¬ 


portant, and may warrant a more detailed analysis. 
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APPENDIX: GOVERNING EQUATIONS 

In this appendix, we sketch the derivation of the Euler- 
Lagrange equations of the energy functional (4), follow¬ 
ing [22]. The variation takes the form 

5£ / \ 

^ = y r-Q [Ns 6Es -b SE^j ds 

-b J ro(^Ms6Ks + ds, (Al) 

where we have introduced the stresses and moments 

Ns = C{Es+ vE^) , N^ = C{nEs + E^) , (A2a) 
Ms =D{Ks + = D{vKs + K^) , (A2b) 

with C = Ehj— and D = Cl'? jVl. (These stresses 
and moments are expressed here relative to the unde¬ 
formed configuration.) 

The deformed shape of the shell is characterised by 
the radial and vertical coordinates r(s) and z(s), as well 
as the angle /3(s) that the normal to the deformed shell 
makes with the vertical direction. These geometric quan¬ 
tities obey the equations [22] 

— = /s cos/3, ^ = /s sin/3, — = jsK,s. (A3) 

We note that one of these is redundant. The varia¬ 
tions SEs,6E^, SKs, 5K^ are purely geometrical, and one 
shows that [22] 


SEs = sec /3 Sr' + fs tan/3 5/3, SE^ = —, 

ro 


SKs = 5/3', 


SK^ = ^5/3. (A4b) 


The variation (Al) thus becomes 


47r 


roNs sec /3 5r -b tqMs 5/3 


— (tqNs sec 13^ - N^jSr ds 

+ J |ro/s3Vstan/3-bM0COS/3-^(^roMs^|5/3ds, (A5) 
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upon integration by parts, whence 

rofsNs tan/3 + cos /3 - ^ (roMs'j = 0, (A6a) 

^ (^roNs sec - N^ = 0. (A6b) 

These equations, together with two of the geometric rela¬ 
tions (A3), describe the shape of the deformed shell. For 
numerical purposes, it is convenient to remove the sin¬ 
gularity at /3 = 7r/2 by introducing the transverse shear 


tension [20, 22], Q = — A^stan/3, expressed here relative 
to the undeformed conhguration. Force balance argu¬ 
ments [20, 22] show that Q obeys 

^ + rofsKsNs + rof^K^N^ = 0. (A7) 

The solution Q = —Ng tan /3 is selected by the boundary 
condition (5(0) = 0. At the poles of the shell, the equa¬ 
tions have singular terms in them, but these singularities 
are either removable or the appropriate boundary values 
are set by symmetry arguments [22]. This allows appro¬ 
priate boundary conditions and values to be derived. 
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